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We briefly summarize the important and critical roles that PHENIX Experiment has 
played in determining the proton’s internal structure in terms of quarks and gluons, and 
their dynamics. Some pioneering measurements by PHENIX Experiment on the motion 
and polarization of quarks and gluons, as well as their correlations inside a fast moving 
proton are presented. Some future opportunities and potentials of PHENIX Experiment 
are also discussed. 
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1. Introduction 

Proton, along with neutron, known as the nucleon, is the fundamental building block of all 
atomic nuclei that make up all visible matter in the universe, including the stars, planets, 
and us. Proton is also known as a strongly interacting, relativistic bound state of quarks and 
gluons (referred as partons) in Quantum Chromo-Dynamics (QCD) [I]. In order to explain 
the origin, the evolution, and the structure of the visible world, it is necessary to explore 
and to understand the proton’s internal structure and its formation. 

The advances of accelerator technology in the last century made it possible to discover the 
quarks, the fundamental constituents of the proton, in the nineteen sixties at Stanford Linear 
Accelerator Center (SLAC) [2, 3], which led to the discovery of QCD. QCD is the theory of 
strong interacting color charges that is responsible for binding colored quarks into the color 
neutral proton by the exchange of gluons [4]. In contrast to the quantum electromagnetism, 
where the force carrying photons are electrically neutral, the force carrying gluons in QCD 
carry the color charge causing them to interact among themselves, which is the defining 
property of QCD responsible for the complex and extremely rich, while mysterious structure 
of the proton, and the countless intriguing phenomena of strong interaction physics. 

Proton structure is not static, but, dynamic, full of features beyond our current knowledge 
and imagination. Quarks and gluons interact strongly when they are far apart - the confine¬ 
ment of color, but, weakly when they are closer - the Asymptotic Freedom; they appear and 
disappear in numbers almost continuously; and they are confined to form the proton while 
moving relativistically [1]. Understanding completely the proton structure and its properties 
(such as its mass, spin, size, and etc.), the emergent phenomena of QCD dynamics, is still 
beyond the capability of the best minds in the world today. It is the great intellectual chal¬ 
lenge to explore and to understand the proton structure and its properties without being 
able to see the quarks and gluons in isolation. 
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Relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory (BNL) was the 
first collider facility in the world to be able to perform experiments with the collisions of 
two relativistic heavy ion beams, as well as two polarized proton beams [5]. From collisions 
of two gold ion beams, RHIC discovered, later confirmed by the experiments at the Large 
Hadron Collider (LHC) at CERN, the strongly interacting quark-gluon plasma (QGP), the 
matter is expected to be only existed in the universe a few microseconds after its birth [6]. 
From collisions of two very energetic polarized proton beams, RHIC, as the first and only 
polarized hadron collider in the world, discovered that from the momentum region accessible 
by RHIC experiments, gluons inside a longitudinally polarized fast moving proton are polar¬ 
ized and have a net positive polarization along the direction of the proton, equal to about 
20% of the proton’s helicity [7, 8]. With one proton beam polarized longitudinally, taking 
the advantage of the parity violation nature of the weak interaction, RHIC has performed 
the best measurements of the net sea quark polarization and its flavor dependence inside a 
polarized proton [9, 10]. With one proton beam polarized transversely, RHIC has provided 
and will continue providing tremendous new opportunities to probe and to study the quan¬ 
tum correlations between the spin direction of the proton and the preference in direction of 
the confined motion of quarks and gluons within the proton [11, 12]. Such quantum corre¬ 
lations, without any doubt, are the most sensitive to the QCD dynamics and the formation 
of its bound states. 

The experimental measurements and discoveries by the PHENIX collaboration, along with 
those by the STAR collaboration and early BRAHMS and PHOBOS collaborations, have 
defined the RHIC science program and will continue to do so. In this short article, we 
briefly summarize the important and critical roles that PHENIX Experiment has played in 
determining the proton’s internal structure, especially emphasizing the recent achievements. 
Some future opportunities and potentials of PHENIX Experiment for exploring the proton 
structure are also discussed. 


2. The helicity structure of the proton 

Parton distribution and correlation functions describe the fascinating relation between a fast 
moving proton in high energy scattering and the quarks and gluons within it. They carry rich 
information on proton’s mysterious partonic structure that cannot be calculated by QCD 
perturbation theory. Parton distribution functions (PDEs) are the simplest of all correlation 
functions, fi/p{x, fj?), defined as the probability distributions to find a quark, an antiquark, 
or a gluon (i = g, q, g) in the proton carrying its momentum fraction between x and x + dx, 
probed at the factorization scale g,. PDEs also play an essential role to connect the measured 
cross sections of colliding proton(s) to the short-distance scattering between quarks and 
gluons. Without them, we would not be able to understand the hard probes, cross sections 
with large momentum transfers, in high energy hadronic collisions, as well as the discovery 
of Higgs particles in proton-proton collisions at the LHC. PDEs are nonperturbative, but, 
universal, and have been traditionally extracted from QCD global analysis of all existing 
high energy scattering data in the framework of QCD factorization[13-16]. The excellent 
agreement between the theory and data on the scaling-violation behavior, or the factorization 
scale /U-dependence, of the PDEs has provided one of the most stringent tests for QCD as 
the theory of strong interaction. 
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Proton has spin 1/2 in the unit of fundamental constant h. That is, proton has a quantized 
finite angular momentum when it is at the rest. Proton is a composite particle made of quarks 
and gluons. Its spin, like its mass and other properties, is an emergent phenomenon of QCD 
dynamics. We do not fully understand QCD if we do not understand how quarks and gluons, 
and their dynamics make up the proton’s spin. 

The proton’s spin 1/2 was once attributed to the superposition of quark spin of three 
valence quarks according to the Quark Model. The proton spin puzzle (also known as the 
spin “crisis”) - a very little of proton’s spin is carried by its quarks, discovered by European 
Muon Collaboration (EMC) about 30 years ago [17, 18], has led to tremendous theoretical 
and experimental activities to explore the proton’s spin structure and to search for an answer 
to one of the most fundamental questions in QCD - how do the spin, the confined motion 
and the dynamics of quarks and gluons inside the proton make up its spin 1/2? 

2.1. The RHIC spin program 

The RHIC spin program, with the only polarized proton-proton collider in the world, was 
designed to probe the internal quark/gluon structure of a polarized proton, and to search 
for the origin of and the answer to the long-standing spin puzzle. The RHIC spin program 
has the following three goals [19]: 

o precision measurement of the polarized gluon distribution g{x, fR) over a large range of 
momentum fraction x to constrain the gluon spin contribution to the proton’s spin, 
o measurements of the polarized quark and anti-quark flavor structure in the proton, and 
o studies of the novel transverse-spin phenomena in QCD. 

Having the collisions of two high energy polarized proton beams with the polarization in 
either longitudinal or transverse direction, the RHIC spin program has a special advantage 
in probing the content of polarized gluons inside a polarized proton and the novel transverse 
spin phenomena in QCD. With its reach in both energy and polarization, the RHIC spin 
program has an unmatched capability to probe the flavor dependence of polarized sea of 
the proton by studying the parity violating single-spin asymmetry of IT^-boson production. 
The RHIC spin program is playing a unique and very critical role in our effort to solve the 
proton spin puzzle. 

Proton is a composite particle of quarks and gluons, and its spin could receive contribution 
from the spin of quarks and gluons, as well as the orbital angular momenta of quarks and 
gluons due to their confined motion. Proton is also a highly dynamic quantum system of 
quarks and gluons full of fluctuations at all distance scales. Numbers of quarks and gluons 
vary at all time, given by the probability distributions to find them, such as the PDFs. 
Similarly, the quark and gluon spin contributions to a polarized proton are given by their 
helicity distributions, 

Afi/p{x, = /Jp(x, - /^^(x, , (1) 

defined as the difference of the probability distribution to find a parton of flavor i with 
the same helicity as the proton, //^, and that of opposite helicity, f~p. The spin 1/2 of a 
longitudinally polarized proton could be decomposed according to the following fundamental 
sum rule [20-23], 

Sp = - = - AE(^^) -I- AG{p^) + [Lq{p?‘) + Lg{p?)) , (2) 
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where AS(^^) and AG{iJ?) are net quark and gluon helicity defined as, 


AS(^^) = f dx A/j/p(x,/U^) = f dx (Art + An + Ad + Ad + As + As) (x,/U^), 

i=q,q-^0 Jo 

AG(^^) = [ dxAfg/p{x,n^)= [ dxAg{x,iJ^). (3) 

do do 


In Eq. (2), the factor 1/2 in the right-hand-side of the equation is the spin of each quark and 
anti-quark; and Lq{ijJ) and Lg{gJ) represent the quark and gluon orbital angular momentum 
contribution, respectively. The parton helicity distributions. A/, are the key ingredients for 
solving the proton spin puzzle. 

The experiment of inclusive electron-proton deep inelastic scattering (DIS), i{k) + h{p) —)■ 
£'{k') -|- X with a large momentum transfer Q = \/—q‘^ = — k'Y 1/fm, discovered 

the quarks inside the proton in the nineteen sixties at SLAC. The same experimental setting 
with longitudinally polarized electron and longitudinally polarized proton can access the 
helicity states of quarks and gluons of the polarized proton. With the approximation of 
one-photon exchange between the scattering lepton and the proton, neglecting the Z-boson 
exchange, the difference of the polarized DIS cross sections with the proton’s spin reversed 
measures the polarized DIS structure function, gi{xB-,Q‘^)-, as 

= (4) 






dx B dx B 


where the terms suppressed by x\{Mp/Q‘^) with proton mass Mp have been neglected, 
xb ^ Q‘^I‘^P • q is the Bjorken variable, and y = q • p/k ■ p is the inelasticity of DIS. In terms 
of QCD factorization [24], the DIS structure function gi can be factorized into a sum of 
parton helicity distributions with perturbatively calculable coefficients expressed in a power 
series of strong coupling ag [25-31]. At the leading order (LO) in ag, 

giixB, = eg [Aq{xB, + Aq{xB, Q^)] , (5) 

where Sg denotes a quark’s electric charge, and the resolution scale of the parton helicity 
distributions is set to be equal to the resolution scale of the exchange virtual photon, p = Q. 
It is clear from Eq. (5) that the structure function gi{xB,Q‘^) is directly sensitive to the 
proton’s spin structure in terms of the combined quark and anti-quark helicity distributions. 
Because of the electroweak probe of the inclusive DIS measurements, the gluon helicity 
distribution Ag enters the expression for gi only at higher order in perturbation theory. Since 
Ag also contributes to the scaling violations (the resolution scale //^-dependence) of quark 
and anti-quark helicity distributions, the access to gluon polarization by DIS measurements 
of gi requires a large level arm in Q^, which could be achieved at a future Electron-Ion 
Collider (EIC) [32]. 

With the various hadronic probes, such as jet(s) and identified hadron(s), and electroweak 
probes, such as the photon, lepton(s), and IT/Z-bosons, the RHIC spin program and its 
measurements of many complementary observables provide a much more direct access to 
gluon polarization, as well as flavor separation of quark polarizations. 
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2.2. The gluon helicity structure 

Longitudinally polarized proton-proton collisions at RHIC allow access to gluon helicity 
distribution, Ag{x,fj?), at LO in perturbative QCD (pQCD). With its detector advantages, 
the PHENIX experiment makes the best connection to Ag via single inclusive tt® production 
at large transverse momentum p± by measuring the inclusive double-helicity asymmetries. 
All, defined as 


All 


Aa 

a 


(^++ - ^+- 
o--r+ + 0-+- ’ 


( 6 ) 


where Aa (a) is the polarized (unpolarized) cross section, and cj++ (ct-i _) represents the 

cross section oi p + p collisions with the same (opposite) proton helicity. In terms of QCD 
factorization, both the polarized and unpolarized cross sections for the single inclusive tt^ pro¬ 
duction can be factorized to show the explicit connection to the parton helicity distributions 
and PDFs, respectively, 


A 


7r» 

LL 


Y.abc^fa/p{xi) ^ 2\fb/p{x2) ^ Aha+b^c+X {xi, X2, Pc) ® Dc^^o{z) 
T,abcfa/p{xi) ® fb/p{x2) ^ aa+b^c+x{xi, X2,Pc) ® Dc^no{z) 


( 7 ) 


where xi and X 2 are momentum fractions of two colliding partons, Dc^n^iz) is the fragmen¬ 
tation function (FF) representing a probability distribution for a parton c of momentum pc 
to fragment into a hadron vr^ of momentum zpc, (8> represents the convolution over the active 
parton’s momentum fractions, and the dependence on the factorization and renormalization 
scales are suppressed. In Eq. (7), Aha+b^c+x (ha+b^c+x) is the polarized (unpolarized) 
scattering cross section of a partonic subprocess, a + b ^ c + X, for two active incoming 
partons of flavor a and b, respectively, to produce a single inclusive parton of flavor c and 
momentum pc. Both polarized and unpolarized partonic cross sections are calculable in 
pQCD, and are available for LO as well as next-to-leading order (NLO) in powers of ag. 
For the kinematic region covered by the RHIC energies, partonic cross sections of gluon 
initiated subprocesses are signihcantly larger than those initiated by quarks and antiquarks. 
That is, measurements of the asymmetry A'^j^ by PHENIX Experiment at RHIC energies 
could provide very sensitive information on the gluon helicity distribution. 

The factorized expression for in Eq. (7) should be valid for all single inclusive hadron 
production so long as the produced hadron transverse momentum is much larger than its 
mass, ph± ^ Mh, to ensure the validity of QCD factorization, and the tt® EEs, Dc.^tto{z), 
are replaced by corresponding hadron EEs, Dc^h{z). This is because all factorized par¬ 
tonic scattering cross sections are insensitive to the details of hadronic states produced. 
Having the universal EEs, Dc^h{z), extracted from other scattering processes, such as 
e"*" -|- e“ —7- h{ph) + X, measurements of A^^ for the production of single inclusive hadron h 
other than provide additional information on the universal parton helicity distributions. 
The PHENIX Experiment measured the asymmetry [8, 33-36], as well as [37]. 

By removing the dependence on the EEs, and replacing the single-parton inclusive cross sec¬ 
tions, Aha+b^c+x and aa+b^c+x by corresponding partonic jet cross sections, 
and , respectively, the factorized expression for in Eq. (7) reduces to the 

factorized expression for which depends on the same parton helicity distributions and 
PDFs. Precise measurements of the double longitudinal spin asymmetry for jet production 
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at RHIC energies could also provide excellent information on parton helicity distributions 
[7, 38-41], 

With the earlier measurements of the double helicity asymmetry, by the PHENIX 

Experiment at the mid-rapidity for y/s = 200 GeV collision energy in 2005 and 2006 [33- 
36, 42], along with measurements of by the STAR Experiment [38-41], the NLO QCD 
global analysis, known as the “DSSV” analysis, concluded that the RHIC data - within their 
uncertainties at that time - did not show any evidence of a net polarization of gluons inside 
the proton [43, 44], 

Recently, based on data collected in 2009 at y/s = 200 GeV at RHIC, which not only 
approximately doubles the statistics of the earlier measurements [35, 37], but also extends 
the range of measured p±, PHENIX Collaboration published new measurements for both 
A'^j^ and Aj^j^ [8], As shown in Eig. 1, three data sets collected during the 2005, 2006, and 



[GeV/c] [GeV/c] 


Fig. 1 (Left) and Aj^j^ (Right) as a function of meson’s transverse momentum pT for 

the 2005 (red circle), 2006 (blue square) and 2009 (black triangle) PHENIX data sets. 

2009 RHIC runs with polarized proton collisions are consistent with each other. 

With the better statistics and extended p± range, new data from 2009 RHIC run could 
naturally put tighter constraints on gluon helicity distribution, Ag{x), and extend the range 
of X over which meaningful constraints can be obtained from QCD global fits. Within exper¬ 
imental uncertainties, the new PHENIX data on and A^j^ are still consistent with the 
existing DSSV analysis, and other recent NLO QCD global analyses of DIS-only data by 
Bliimlein and Bottcher (BBIO) [45] and Ball et. al. (NNPDF) [46, 47], as well as the analysis 
by Leader et. al. (LSSIO) [48] based on both DIS and Semi-Inclusive DIS (SIDIS) data. Since 
various analyses used different assumptions on the symmetry properties of parton helicity 
distributions, and different functional forms for the distributions at the input scale, the 
determination of the gluon helicity distribution, Ag{x, p?) varies, see the detailed discussion 
in the PHENIX publication [8]. 

Although the new PHENIX data on AJ*^ and A^^ do not show any significant asymme¬ 
try [8], the new STAR data on A'^^^ of inclusive jet production from the same 2009 run at 
RHIC show a non-vanished double-spin asymmetry over the whole range 5 < pr ^ 30 GeV 
[7], which differs from the previous results on jet production published by the STAR Collab¬ 
oration [38-40]. The most resent DSSV QCD global analysis, referred as “DSSV-i-i-” [49], 
shows that the new PHENIX data and STAR data are actually consistent with each other. 
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as shown in Fig. 2(Left). The new data sets and new DSSV++ QCD global analysis lead to 
a significant net positive gluon helicity contribution from the x region probed by the current 
measurements at RHIC, 0.05 < x < 0.2, with a better constrained profile for the global 
analysis, as shown in Fig. 2(Right). For the first time, we have a clear experimental evidence 
to show that gluon could give a non-vanish contribution to the proton’s spin, at least from 
the momentum region probed by the current RHIC measurements. 


ji° (GeV/c) 




Fig. 2 Left: Double-spin asymmetry All of new PHENIX and STAR data are compared 
with the new DSSV+-(- QCD global analysis [49]. Right: the prohle for the integrated 
gluon helicity contribution in the x region currently probed at RHIC for ^/s = 200 GeV. 

With new measurements of vr® in large forward rapidity region at = 500 GeV, a higher 
center-of-mass energy of polarized proton-proton collisions at RHIC, the PHENIX Experi¬ 
ment with its forward calorimeter can access the gluon helicity distribution at a momentum 
fraction x as small as 0.002, which will considerably improve our knowledge of the gluon 
spin contribution to the proton’s spin in the x-range of 0.002 — 0.2 [19]. 

2.3. The sea quark helicity structure 

The combination of the large body of inclusive DIS data has provided excellent measurements 
of combined quark and antiquark helicity structure, and established that the up quarks and 
anti-upquarks combine to have net polarization along the proton spin, whereas the down 
quarks and anti-downquarks combine to carry negative polarization. The “total” Au + Au 
and Ad + Ad helicity distributions are very well constrained for the medium to large x, 
contributing to about 30% of the proton’s spin. 

However, the light sea quark and anti-quark structure is still far from being well-understood 
[50], and its helicity distributions still carry large uncertainties, even though there are some 
constraints by the SIDIS data. Taking the advantage of the fact that W bosons couple 
only the left-handed quarks and right-handed antiquarks {uLdn —>• W~^ and dLUR —)• W~), 
measuring the single longitudinal spin asymmetry, Al = (ct+ — a-)/{a+ -|- fT_), of the parity 
violating production of lE-bosons from flipping the helicity of one of the polarized colliding 
protons at RHIG can probe the flavor dependence of sea quark helicity distributions, Aq 
and Aq. Eor example, the difference of the W~^ yield from the diagram on the left in Eig. 3 
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Fig. 3 Sketch for W~^ production from a leading order annihilation of a right-handed 
anti-downquark and a left-handed up quark in the collision of a polarized proton (top) and a 
unpolarized proton (bottom). The produced W~^ then decays to a positron and an electron 
neutrino. 

and that from the diagram on the right is of good information on anti-downquark helicity 
distribution, A(i(xi) = d'^{xi) — d^(xi), where the subscript -|- (—) indicates the positive 
(negative) helicity state of the polarized proton, the superscript -|- represents the handiness 
of the anti-downquark, and xi is the momentum fraction of the polarized proton carried 
by the anti-downquark. Production of the VP-boson requires a large momentum transfer 
at a scale where higher order QCD corrections can be evaluated reliably, and is free from 
uncertainties in fragmentation functions that are needed for probing light hadrons. The 
measured VP-boson cross sections in unpolarized hadronic collisions at Tevatron and the 
LHC, as well as at RHIC, have confirmed our theoretical understanding of the production 
mechanism. 

The measurements of of VP^-bosons at RHIC as a function of decay lepton’s transverse 
momentum and rapidity provide the uniquely clean information on the sea quark helicity 
distributions, as well as its flavor structure inside a polarized fast moving proton. With the W 
mass setting up the hard scale for the momentum transfer of the collisions, RHIC’s IP-physics 
program could probe the quark and anti-quark helicity distributions at a medium momentum 
fraction, 0.05 < x < 0.4, providing strong constraints on the sea quark contribution to the 
proton’s spin. Even more important, RHIC’s IP-physics program could provide the unique 
and much needed information on the sea quark helicity asymmetry, such as Ad{x) — Au(x), 
in a perfect kinematic region, in view of the puzzling sign change of unpolarized sea quark 
distributions, d{x) — u{x), at x ~ 0.3 observed by the E886/NuSea Collaboration at Fermilab 
in its extraction of d{x) — u{x) over the region 0.02 < x < 0.345 from the measurements 
of Drell-Yan process [51-53]. The sea quark asymmetry in both momentum and helicity 
distributions is fundamentally important for understanding the proton structure. 

From the background subtracted yields of positions and electrons in p-p collisions at ^/s = 
500 GeV at RHIC, PHENIX experiment extracted the IP^ and IP“ boson production cross 
sections, a{pp^W~^X) x BR{W~^ ^e'^Ve) = 144.1 ± 21.2(stat)15Q'^3(syst) ± 15%(norm)pb, 
and a{pp^W~ X) x BR{W~ ^e~i7e) = 31.7 ± 12.1(stat)lg°2^(syst) ± 15%(norm)pb, respec¬ 
tively, where BR is the IP’s decay branching ratio [10]. The measurements of PHENIX 
Collaboration, together with that of STAR Collaboration [54], provided in fact the very 
first IP-boson production cross sections in p-p collisions, which are consistent with earlier 
measurements in p-p collisions at CERN and Fermilab [10]. 
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Fig. 4 PHENIX measurement ol Al for electrons (Left) and positions (Middle) from W 
and Z decays [10], along with theory predictions [55]. Right: Expected uncertainties for Ai 
of production for PHENIX and STAR after the 2013 run. The asymmetries have been 
randomized around the central value of DSSV analysis. 

With its early data on leptonic W decay in polarized p-p collisions at y/s = 500 GeV, 
RHIC’s IT-physics program started to impact our understanding of proton’s polarized sea 
structure [9, 56]. In Fig. 4, the very first PHENIX measurement of the parity violating single 
spin asymmetry, A^, of electrons (Left) and positions (Right) from W and Z decays are 
presented [10]. Although the experimental uncertainties are clearly large, the asymmetries 
are consistent with the early theory expectations. 

With the new STAR data on A^^ [9], and the expected experimental uncertainties for 
Ai measurements of both PHENIX and STAR from the recent polarized p-p runs at RHIC 
[19], as shown in Fig. 4 (Right), new and updated DSSV global analyses of parton helicity 
distributions clearly indicate that the sea quark structure shows hints of not being SU(2)- 
flavor symmetric: the Au distribution has a tendency to be mainly positive, and the Ad 
anti-quarks likely carry an opposite polarization. From the profiles of the global fits in 



Fig. 5 Expected profiles from DSSV global analyses of parton helicity distributions 
for Am (left) and Ad (right): DSSV-t with new SIDIS data, and DSSV-t-1- with new STAR 
data or with projected data on of boson production in Fig. 4 (Right). 
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Fig. 5, the new RHIC data on the parity violating single spin asymmetries will provide a 
powerful constraint on the sea quark structure and its contribution to the proton’s spin. 

QCD global analyses of all available data from DIS, SIDIS, and RHIC have provided 
good information on quark/gluon helicity structure of a polarized fast moving proton, and 
quark/gluon spin contribution to the proton’s spin [32]. From the region of momentum 
fractions accessible by existing experimental measurements, the net helicity of quarks and 
antiquarks combined could account for about 30% of proton’s spin, while gluons are likely 
to give a positive contribution to the proton’s spin, close to 20% from the most recent RHIC 
data, but with a hugh uncertainty from the region of gluon momentum fractions that is 
not available to all existing and near future facilities [49]. While we need the hnal word of 
the gluon spin contribution to the proton’s spin from the proposed future EIC [32], it is 
fundamentally important to investigate the conhned transverse motion of quarks and gluons 
inside a fast moving proton, and their contribution to the proton’s spin. 

3. Transverse parton structure of the proton 

To have a direct access to the transverse parton structure of a fast moving proton in high 
energy scattering requires a new type of physical observables - the controlled probes that 
have at least two distinctive momentum scales: a large scale Q ^ Aqcd ~ 1/fm to localize 
the probes to “see” the particle nature of quarks and gluons, and a small scale q Q to 
be sensitive to the conhned transverse motion or spatial distribution of quarks and gluons 
inside the proton. 

Theoretically, the simplest transverse parton structure of the proton could be characterized 
by the Wigner distributions Wiip{x, kx, bx) - the probability distributions to hnd a parton 
of havor i{= q, q,g) inside the fast moving proton, carrying its momentum fraction between 
X and X + dx and having a transverse momentum kx at a transverse spatial position bx 
with respect to its center [57], More complex transverse parton structure of the proton could 
be represented by multi-parton correlation functions in terms of longitudinal momentum 
fractions, transverse momenta and positions of all active partons. Although the Wigner 
distributions are used extensively in other branches of physics [58], there has not been any 
known way to measure the quark and gluon Wigner distributions in high energy experiments. 

However, various reductions of the Wigner distributions could be probed in high energy 
experiments [32]. By integrating over bx, the Wigner distributions are reduced to the three- 
dimensional (3D) momentum distributions of quarks and gluons inside the fast moving 
proton - transverse momentum dependent PDFs (TMDs), fi/p{x, kx)', and the TMDs could 
be further reduced to the better known PDFs, fi/p{x,g?) by integrating over kx with a 
proper ultraviolet (UV) renormalization when kx —?• oo for the quark/gluon operators defin¬ 
ing the PDFs at the renormalization scale fi. Similarly, by integrating over kx, the Wigner 
distributions are reduced to the impact parameter distributions, fi/p{x,bx), the probability 
distributions to hnd a parton of havor i and momentum fraction x at a transverse spatial 
position bx, which are effectively the tomographic images of quarks and gluons inside the 
proton. Furthermore, by taking the Fourier transform of bx into the conjugate momentum 
A.T, fi/p{x, bx) could be transformed into the generalized PDFs (GPDs), such as Hi/p{x, t) 
and Ei/p{x,^,t) at ^ = 0 with t = — [59]. The GPDs could be extracted from measure¬ 

ments of exclusive processes in lepton-hadron DIS or in ultra-peripheral hadronic collisions. 
The TMDs and GPDs represent various aspects of the same proton’s transverse structure of 
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quarks and gluons that could be probed in high energy scattering [32]. In fact, some GPDs 
are intimately connected with the orbital angular momentum carried by quarks and gluons 
[60]. The Ji’s sum rule is one of the examples that quantify this connection [21], 

1 

'^1=0 M Wq{x, t) + Eg(x, C, t)] , (8) 

Z t^U Jq 

which represents the total angular momentum Jq (including both helicity and orbital con¬ 
tributions) carried by quarks and anti-quarks of flavor q. A similar relation holds for 
gluons. 

With the polarized high energy proton beam(s), the RHIC spin program could access the 
proton’s transverse structure of quarks and gluons by measuring various TMDs, as well 
as GPDs via exclusive processes in ultra-peripheral p-p and p-A collisions. The RHIG spin 
program could measure many emergent QCD phenomena, such as quantum correlations 
between the proton’s spin direction and preference of the confined transverse motion of 
quarks and gluons inside the polarized proton, known as the Sivers effect [61]. Furthermore, 
it could probe quantum correlation between the spin direction of the fragmenting parton 
and preference in direction where the produced hadron emerges, known as the Gollins effect 
[62]. 

3.1. Transverse single-spin asymmetry - A 

Transverse single-spin asymmetry (SSA), A^r = (cj(st) — <7(—ST))/(o'('Sr) + c’'(—sr)), is 
defined as the ratio of the difference and the sum of the cross sections when the spin of 
one of the identified hadron st is flipped. Large SSAs of inclusive single pion production 
with a large momentum transfer in hadronic collisions have been consistently observed since 
nineteen seventies, as shown in Fig. 6, although the asymmetries were once thought impos¬ 
sible in QGD [63]. With over two decades of intense theoretical as well as experimental 
efforts, our understanding of the observed large SSAs in QCD has been much improved. 
Large SSAs are not only possible in QCD, but also carry extremely valuable information 
on the transverse motion and structure of quarks and gluons inside a transversely polarized 
proton. 



Fig. 6 Persistent transverse single-spin asymmetry, A^r, in hadronic single inclusive pion 
production. 

From the parity and time-reversal invariance of QCD dynamics, non-vanish transverse 
SSAs are necessarily connected to the transverse momentum of quarks and gluons inside the 
transversely polarized proton [64]. Two complementary QCD-based approaches have been 
proposed to describe the physics behind the measured SSAs: the TMD factorization approach 
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[61, 62, 65, 66] and the collinear factorization approach [67-72]. In the TMD factorization 
approach, the asymmetry is attributed to direct correlation between spin and transverse 
momentum, which are represented by the TMD PDFs or FFs, such as the Sivers and Collins 
functions, respectively. On the other hand, the asymmetry in the collinear factorization 
approach is generated by twist-3 collinear PDFs or FFs. With the transverse momenta 
of all active partons integrated, the explicit spin-transverse momentum correlation in the 
TMD approach is now replaced by the integrated net effect of spin-transverse momentum 
correlations, which are effectively generated by QCD color Lorentz force [73]. 

Two approaches each have their own kinematic domain of validity, while they are con¬ 
sistent with each other in the perturbative regime where they both apply [74-78]. The 
TMD factorization approach is more suitable for evaluating the SSAs of scattering processes 
with two very different momentum transfers, Q ^ q ^ AqcDj while the collinear factoriza¬ 
tion approach is more relevant to the SSAs of scattering cross sections with all observed 
momentum transfers hard and comparable: Q' ^ Q ^ Aqcd- In hadornic collisions, the sin¬ 
gle inclusive hadron production at high is better treated in the collinear factorization 
approach, while the Drell-Yan and WjZ production at a low transverse momentum, <C Q, 
needs the TMD factorization approach. 

3.2,. An of single hadron production 

The inclusive single hadron production at large transverse momentum pT at RHIC has 
effectively one large momentum transfer at 0{pt) Aqcd- Collinear factorization approach 
is more suited for studying the transverse SSAs of the inclusive single hadron production. 
In terms of the collinear factorization, the transverse SSAs are effectively power suppressed 
observables comparing to the production cross sections. That is, the asymmetries are in 
general small, except in the region of phase space where the momentum spectrum of the 
observed pions is very steep, and a small shift in the spectrum could make the difference 
of two cross sections with the proton spin flipped to be comparable with the cross section 
itself, leading to a significant value of the asymmetry. 

2 0.2 
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Fig. 7 as a function of xp at ^/s = 62.4 GeV (Left) measured by PHENIX Experiment, 

and at various center-of-mass energies in hadronic collisions (Right). 

The left plot in Fig. 7 shows at ^/s = 62.4 GeV measured by PHENIX Collaboration 
with the PHENIX detector at RHIC [12]. While there is a significant, nonzero asymmetry 
rising with xj? > 0 in the forward direction of the polarized proton beam, no such behavior 


P+P ^ 71 " + X 

* PHENIX 7i“3.1<ti<3.8S=62.4 GeV 
A E704 7t“ (s=19.4GeV 

* STAR 7i“ <t|>= 3.3, fc200 GeV 
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can be seen at negative xp < 0 where the asymmetries are consistent with zero. The typical 
transverse momentum of these data is relatively low, over the range of 0.5 < pp < ^ GeV 
[12]. In Fig. 7 (Right), the PHENIX data of at ^/s = 62.4 GeV are compared with the 
measurements by Fermilab E704 at ^/s = 19.4 GeV and those by STAR Collaboration at 
^/s = 200 GeV. Although these measurements were carried out with slightly different detec¬ 
tor acceptances, there is a general agreement between the xp dependence of nonvanishing 
asymmetries. The observed asymmetries appear to be independent of the center-of-mass 
energy of the collisions, which varies over one order of magnitude from ^/s = 19.4 GeV to 
200 GeV. The data set plotted in Eig. 7 (Right) also covered an interesting range of trans¬ 
verse momentum of the observed pions from a semi-hard scale at pp ~ 0.5 GeV to a hard 
scale at Pp > 2 GeV for the STAR data [11]. 



Fig. 8 as a function of the cluster’s Feynman parameter xp (Left) and transverse 

momentum pp (Right) measured by PHENIX Collaboration at ^/s = 200 GeV at RHIC [12]. 


By resolving the two electromagnetic clusters from the two photons of vr^ —)• 7 -|- 7 decay, 
PHENIX detector is capable of detecting with its energy, < 20 GeV. However, with 
increasing energy, the opening angle between the two photons becomes so small that their 
electromagnetic clusters fully merge in the PHENIX detector. The reconstruction of from 
the two-gamma decay mode limits the xp range to below 0.2 at ^/s = 200 GeV, where 
transverse SSAs are small. To overcome this limitation, the data analysis of transverse SSAs 
is also done for inclusive electromagnetic clusters [12]. The clusters are dominated from 
7 + 7 ) with some contributions from direct and other photons and decay of 7 and 
other charged particles. The detailed discussion on the decomposition of the clusters can be 
found in Ref. [12]. 

The left plot in Fig. 8 summarizes the XiT’-dependence of the cluster Apf for two differ¬ 
ent pseudorapidity ranges. Within statistical uncertainties the asymmetries in the forward 
direction of the polarized proton Apf rises almost linearly with xp, while the asymmetries in 
the backward direction xj? < 0 are found to be consistent with zero. The non-vanish asym¬ 
metries in the forward region are of similar size compared to the vr^ asymmetries at other 
center-of-mass energies as shown in Fig. 7 (Right), which should not be too surprising since 
the inclusive electromagnetic clusters in this momentum regime are mainly formed by the 
photons from vr^ decay. 
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The right plot in Fig. 8 presents the cluster yljV; as a function of transverse momentum 
Pt for values of \xf\ > 0.4 [12]. The asymmetry rises smoothly and then seems to saturate 
above px > 5 GeV/c. When px = 0, Ajy should vanish from the symmetry, and is expected 
to rise with px when px is relatively small. Since transverse SSAs of the single inclusive 
hadron production are intrinsically a power suppressed observable at large px, the Aj\f is 
expected to decrease when px becomes sufficiently large [69, 74]. Explicit calculations of A at 
in collinear factorization approach indicates, A^f oc px/u, and consequently, the asymmetry 
does not fall as quickly as 1/px in the forward region when xp > xx = ‘Ipxj^fs [71]. The 
Px dependence of the Api should be sensitive to the underline dynamics generating the 
transverse SSAs. Again, asymmetries in the negative xp region are found to be consistent 
with zero within statistical uncertainties. The transverse SSAs of single inclusive hadron or 
inclusive cluster measured by PHENIX Experiment are consistent with those measured by 
STAR collobration at RHIC. 


3.3. An of Drell-Yan (j* and WfZ) production 

The observed transverse SSAs of single inclusive hadron production could be generated by the 
spin-motion correlations and/or coherent multiple interactions before and/or after the hard 
collision [71]. That is, in the terminology of TMD factorization approach, both the Sivers- 
like (before) and the Collins-like (after) effects are responsible for the observed asymmetries. 
Assuming the correlation before the hard collision dominates the observed asymmetry in the 
inclusive pion production could actually lead to a sign puzzle for the relation between the 
Sivers function and the corresponding twist-3 quark-gluon correlation functions [79]. Appar¬ 
ently, the Collins-like correlations after the hard collision is very important for understanding 
the transverse SSAs in p-p collisions, as well as resolving the sign puzzle [80]. Without being 
able to disentangle different effects, it is difficult, if not impossible, to fully understand the 
physics and the mechanism to generate the novel phenomena of transverse SSAs. 

Unlike the single inclusive hadron production in hadronic collision, Drell-Yan massive 
lepton pair production, either via a virtual photon or a PF/Z boson, could have two natural 
and distinctive momentum scales: the invariant mass of the lepton pair, Q (~ My/jz) and 
the transverse momentum of the pair, qp. The most events have the strong ordering of these 
two scales, Q » qp., which is natrual for applying the TMD factorization approach to probe 
the transverse motion of quarks directly. Since the vector bosons or their decay leptons 
do not interact strongly after they are produced, Drell-Yan process is ideal for study the 
initial-state TMDs, such as Sivers functions [50]. 

One of the key differences between the QCD collinear factorization approach and QCD 
TMD factorization approach to the high energy scattering process is the universality of 
the non-perturbative long-distance physics. Within the collinear factorization approach, all 
non-perturbative long-distance physics are represented by universal (process independent) 
hadronic matrix elements whose operators are made of quark/gluon correlators that are effec¬ 
tively localized in space to the size of the hard collision i/Q <C l/fm. The non-perturbative 
long-distance physics in the TMD factorization approach is represented by TMDs. The exist¬ 
ing definitions of TMDs involve quark-gluon correlators with non-local gauge links covering 
infinite size in space [81]. It is the non-local nature of the TMDs that led to their potential 
process-dependence, such as the sign change of the Sivers function extracted from SIDIS 
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in comparison with that extracted from the Drell-Yan process [66]. This non-universality 
is a fundamental prediction from the gauge invariance of QCD and the TMD factorization 
approach. The experimental test of this sign change is one of the open questions in hadronic 
physics and will provide a direct verification of QCD TMD factorization. 





Fig. 9 Ajsf of Drell-Yan like process with a reconstructed W~^ (Left), W~ (Middle), or a 
virtual photon (Right), as a function of vector boson’s rapidity in p-p collisions at RHIC. 
Also plotted are theoretical predictions based on Sivers functions extracted from low energy 
SIDIS measurements [82]. 


With the high energy polarized proton beams, RHIC can produce W and with rela¬ 
tively low transverse momentum qt The transverse spin program at RHIC has 

the unique advantage to measure the Sivers effect in the Drell-Yan like process at a true 
hard scale - mass ofW/Z bosons. Figure 9 shows the expected uncertainties for transverse 
SSAs of reconstructed (Left, Middle) and Drell-Yan production (Right) from STAR and 
PHENIX measurements [19]. Also plotted curves are theoretical predictions based on our 
knowledges of Sivers functions from SIDIS measurements [82]. These asymmetries provide 
an essential test for the fundamental QCD prediction of a sign change of the Sivers function 
in hadronic collisions with respect to that in SIDIS, as well as a test of our understanding 
of scale dependence (or evolution) of the Sivers functions. 

Transverse SSAs and their xp and pp dependence in hadronic collisions are well- 
established, and they are emergent phenomena QCD dynamics. Along with the transverse 
spin programs at JLabl2 and at COMPASS, the transverse spin program at RHIC with 
its PHENIX and STAR detectors opens up a new domain of QCD dynamics, sensitive to 
the confined motion, as well as the quantum correlation between the motion and intrinsic 
particle property, such as spin, of quarks and gluons inside a QCD bound state, like the 
proton. It also provides controllable and promising access to the transverse structure of the 
proton. 

4. Summary and outlook 

Proton is a dynamic system of confined quarks and gluons. Proton itself is an emergent 
phenomenon of QCD dynamics. Proton structure is dynamic as well. Probability distribu¬ 
tions to find quarks and gluons, and their correlations inside the proton are fundamental. 
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corresponding to various aspects of the dynamic and complex system of the proton. Our 
knowledge of proton’s internal structure is still very limited, although we have learned a 
lot with the advances of accelerator and detector technologies, as well as theoretical break¬ 
throughs to identify the controllable probes to “connect” the quarks and gluons inside the 
proton to the leptons and hadrons observed in modern detectors. 

RHIC, with the only polarized proton-proton collider in the world and the PHENIX and 
STAR detectors, has helped address open questions concerning the proton structure, and 
identified new open questions, puzzles and challenges. The PHENIX Experiment, along with 
STAR Experiment, has played critical roles in defining the RHIC physics program and in 
extracting valuable information for determining proton’s internal structure. The PHENIX 
Experiment helps determine the helicity structure of the proton, in particular, the extrac¬ 
tion of gluon and sea quark helicity structure inside a polarized fast moving proton, by 
measuring An of single inclusive hadron production, as well as the parity violating Al oi 
production. By measuring Atv of many observables, including the production of 7r°, 
rj, direct photon, WjZ bosons, J/V', and the innovative inclusive electromagnetic clusters, 
the PHENIX Experiment help define and explore the richness and excitements of the RHIC 
transverse spin program. 

With the upcoming p-A run at RHIC, the PHENIX Experiment will help open a new 
frontier of QCD to explore the fundamental nuclear structure in terms of quarks and gluons. 
With the help of the polarized proton beam, the PHENIX Experiment could search for the 
imprint of spin correlations among the nucleons inside a nucleus, without seeing nucleon, 
but, quarks and gluons. With the future upgrades of the PHENIX detector to sPHENIX 
and to fsPHENIX, the PHENIX Experiment could be even more powerful in delivering new 
and accurate measurements helping explore the internal structure of the proton. 
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